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Standing waves near the zigzag and armchair edges, and their Berry's phases are investigated. 
It is suggested that the Berry's phase for the standing wave near the zigzag edge is trivial, while 
that near the armchair edge is non-trivial. A non-trivial Berry's phase implies the presence of a 
singularity in parameter space. We have confirmed that the Dirac singularity is absent (present) 
in the parameter space for the standing wave near the zigzag (armchair) edge. The absence of 
the Dirac singularity has a direct consequence in the local density of states near the zigzag edge. 
The transport properties of graphene nanoribbons observed by recent numerical simulations and 
experiments are discussed from the point of view of the Berry's phases for the standing waves. 



The electronic property of a graphene nanoribbon dif- 
fers greatly from that of a carbon nanotube. A metal- 
lic carbon nanotube exhibits a high mobility, while a 
graphene nanoribbon shows a transport gap. [lHJ] The 
high mobility observed in metallic carbon nanotubes in- 
dicates that the scatterers are not effective in produc- 
ing backward scattering. There should be a mechanism 
which suppresses the backward scattering in metallic car- 
bon nanotubes. 

Suppose that an electron with momentum k is coming 
into the impurities which are represented by the circles 
in Fig. []Ja). The electron is scattered by the impurities 
and changes its momentum direction. Let us consider 
the probability amplitude that the electron is scattered 
in the backward direction, as shown in Fig. HJa) . In this 
case, the final wave function, $_k, is given by rotating 
the wave vector of the initial state, <&k, by — it, so that we 
have the relationship between the initial state and final 
state as 



$_ k = U(— 7r)$ k , 



(1) 



where U(6) is a rotational operator with angle 6. Note 
that this particular path shown in Fig. HJa) is not the 
unique path that an electron can follow. There is an 
another path that an electron can follow, which we denote 
it by the lines in Fig. QJb). The new path relates to the 
original path through the "time reversal" . We denote the 
final state in this "time reversal" path by <&'. Because 
the final state is given by rotating the wave vector by 
+7r, we have the relationship, 



$'_ k = C/(+7T)$ k , 



(2) 



between the initial state and the final state. By eliminat- 
ing the wave function of the initial state from Eqs. (p} 



and we get 



$'_ k = £/(2tt)$- 



(3) 



Now, the total backward scattering amplitude is given 
by the sum of $_ k and $'_ k as <J> bs = $'_k + *-k = 
[1 + U(2tt)] $_k- Because the wave function gets an ex- 
tra phase shift of — tt (called the Berry's phase [H, Hj]) 
through a rotation of the wave vector around the Dirac 
point, that is, U(2ir) is equivalent to —1 (U(2tt) = e~ l7r ), 
a "time-reversal" pair of backward scattered waves can- 
cels with each other, i.e., $bs = Q. 



(a) 



(b) 



FIG. 1: The finial state of the scattering path shown in (a), 
"l?_k, combines with the final state of the scattering path 
shown in (b), <J?'_ k , to suppress the total backscattering am- 
plitude, <I>_k + ^'-ki due to the Berry's phase. The two scat- 
tering paths, (a) <&k —¥ ■ ■ ■ — > ^-k and (b) <E> k — > • • • — > *'-k; 
are related with each other by "time-reversal". 



This, absence of backward scattering mechanism pro- 
vides us with a simple solution explaining the high mobil- 
ity, and the existence of a singularity at the Dirac point 
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is essential to the nontrivial phase shift of tt. However, 
it is not obvious whether the wave functions realizing in 
graphene nanoribbons can acquire a nontrivial Berry's 
phase or not. The eigen state near the edge is the stand- 
ing wave resulting from the interference between an in- 
cident wave and the edge reflected wave. If the Berry's 
phase for the standing wave is trivial, then the scatters, 
such as a potential disorder created by charge impurities, 
may give rise to backward scattering. As a result, the mo- 
bility of a graphene nanoribbon decreases considerably 
than that of a carbon nanotube. [3-11] In the present pa- 
per, we study the Berry's phases of the standing waves 
near the zigzag and armchair edges, and their effects on 
the local density of states and the transport properties 
of graphene nanoribbons. 

In Fig. [21 we consider the zigzag edge parallel to the 
x-axis, by which translational symmetry along the y-axis 
is broken. Thus, the incident state with wave vector 
[k x , k y ) is elastically scattered by the zigzag edge, and the 
wave vector of the reflected state becomes (k Xl —k y ). By 
contrast, the armchair edge parallel to the y-axis breaks 
translational symmetry along the cc-axis, so that the wave 
vector of the reflected state is (— k x ,k y ). Since the Bril- 
louin zone (BZ) is given by rotating the hexagonal lattice 
by 90°, one can see in Fig. [2] that for the incident state 
near the K point, the reflected state by the zigzag edge 
is also near the K point, while the reflected state by the 
armchair edge is near the K' point. Hence, the scattering 
by the zigzag edge is intravalley scattering, while that by 
the armchair edge is intervalley scattering. 



y 




FIG. 2: A hexagonal unit cell of graphene consists of A (solid 
circle) and B (open circle) atoms. The xy coordinate system 
is fixed as shown, (right) The hexagon represents the first 
BZ of graphene, and the corners of the hexagon are K and K' 
points. 

To begin with, let us consider the scattering problem 
for the zigzag edge. Since the zigzag edge is not the 
source of intervalley scattering, we focus on only the elec- 
trons near the K point. The incident and reflected waves 
are represented by the Bloch states. The Bloch state in 
the conduction energy band is written as 

e »k-r i / -i0(k)\ 

S£ )k (r) = ^ c K ,k, <S>km = ^\ i )> W 



where V is the area of sample, k is the momentum mea- 
sured from the K point, and 9(k) is the polar angle be- 
tween the vector k and the fc^-axis. Note that the Bloch 
function acquires a nontrivial Berry's phase of — tt by a 
rotation of the wave vector around the K point as 

fj ^KM\~WK*)de = -TT. (5) 

The standing wave modes arise from the interference 
between the incident waves and the reflected waves. Let 
the wave vector of the incident wave be k — (k x ,k y ). 
Then the vector of the elastically reflected wave is given 
by k' = (k x ,—k y ), and the standing wave near the zigzag 
edge is written as 

^,k(r) = ^, k (r) + e-^ ik ,(r), (6) 

where the phase e IZ should be determined in such a man- 
ner that *L k (r) [= *(*K,k,A( r )'*K,k,B( r ))] satisfies the 
boundary condition [8| . We take the following boundary 
condition for the zigzag edge 0, 

I™ *K,k, B (r) = 0, (7) 

which describes a situation in which the B-atoms located 
slightly above the horizontal dashed line at y = in Fig. [2] 
are separated from the lower semi-infinite graphene for 
y < 0. The boundary condition of Eq. (J?) is satisfied 
when e lz — —1, and the standing wave is written as 

where L x is the length of the zigzag edge, and N(y) = 
N for y < and N(y) — otherwise. The value 
of N is determined by the normalization condition. We 
note that the standing wave of Eq. ([6]) reproduces the 
result of tight-binding lattice model (see Appendix B of 
RefE3). 

First, we examine the Berry's phase of the standing 
wave. For simplicity, let us consider the case 0(k) = 
7r/2 for Eq. ©. This corresponds to a normal incident 
process to the zigzag edge [k x — and k y > 0]. The 
amplitude for this process does not vanish as 

Hie M - N (yU- C0S ( k vV)\ (q\ 
* K - k * {r) -7m\sin(k y y) J' (9) 

This standing wave includes a backward scattering ampli- 
tude, and the existence of the backward scattering am- 
plitude indicates that the Berry' phase of the standing 
wave vanishes. Indeed, since 0(k') = — 6(k), the Berry's 
phase of the reflected state is given by +ir as 

/ 2T (fcl-|rlfc>^ = +^ (10) 



3 



while that of the incident state is — tt [see Eq. (J^J]. The 
Berry's phase of — tt for the incident state is canceled 
by the Berry's phase of tt for the reflected state, and 
hence the Berry's phase for the standing wave vanishes in 
total. The trivial Berry's phase can be understood from 
Eq. ©, in which the Bloch function of the standing wave 
is real [6|. 

Next, we study the standing wave near the armchair 
edge. The K and K' points need to be considered simul- 
taneously in the case of armchair edge since the armchair 
edge is the source of intervalley scattering. Suppose that 
the wave vector of the incident wave measured from the 
K point is k = (k x ,k y ). Thus, the wave vector of the 
elastically reflected wave measured from the K' point is 
given by k' = (— k x , k y ), and the standing wave is written 
as 



*k(r) = *K. k (r) + e m <i>^, k ,(r). 



(11) 



It can be shown that the Bloch state near the K' point 
is given by $K',k'( r ) = (e^'V^^K.k & Thus, we 
may write 



*k(r) 



e ik y y / e +ik^x 



(12) 



where N(x) = N ^ for x < and iV(ic) = other- 
wise. The relative phase e la can be determined by the 
boundary condition for armchair edge. The value of e ta 
must be —i, which will be derived elsewhere. [l2| Note 
that the Bloch functions for the K and K' points are the 
same. Thus, the Berry's phase of the standing wave near 
the armchair edge is given by — tt. As a result, we can 
expect that the absence of backward scattering mecha- 
nism is valid near the armchair edge. Note, however, 
that the absence of backward scattering mechanism for 
the armchair edge is not identical to that discussed so 
far for carbon nanotubes. This is because the notion of 
a "time-reversal pair" of scattered waves in the original 
argument [7| should be replaced to a true time-reversal 
pair of scattered waves near the armchair edge. The time- 
reversal state of Eq. (|T2|) is given by 13 1 



-ik y y 



-ia +ik x x 



(13) 



A magnetic field breaks time-reversal symmetry, so that 
it invalidates the new absence of backward scattering 
mechanism and leads to weak anti-localization. It is also 
interesting to point out that the Berry's phase for the 
standing wave near the armchair edge is robust against 
an ordinary mass term, while it is not robust against a 
topological mass term [141 ] . 

The appearance of a nontrivial Berry's phase is re- 
lated to a singularity in parameter space 0, |f|. Here, 
we show the relationship between the singularity at the 
Dirac point (k = 0) in the momentum space and the 



standing wave. For the zigzag edge, k x is a continuous 
variable, and therefore the spectrum with k y = may 
cross the Dirac singularity. However, Eq. (J5J shows that 
the standing wave with k y = does not exist, that is 
k =o( r ) = 0- Thus, the Dirac singularity seems to be 
separated from the spectrum of the standing wave. To 
explore this point further, we consider a zigzag nanorib- 
bon by introducing another zigzag edge at y = —L, in 
addition to the zigzag edge at y = 0. Suppose that the 
edge atoms at y = —L arc B-atoms, which imposes the 
boundary condition on the wave function at y = —L 



as lim,. 



equation for (k x ,k y ), 



0. This leads to the constraint 



k y L + 0(k) = 7T.7T, 



(14) 



where n is an integer. We note that this equation pro- 
vides the solutions of —k x = k y / tan(fcy L) which was 
obtained by Brey and Fertig in Ref. [the minus sign in 
front of k x is a matter of notation] . Note that n must be 
a nonzero integer because the equation does not possess 
a solution when n = 0. In Fig. [3l we plot the solutions 
of Eq. (|14|) for the cases n = 1 and 2. The solution for 
n > 2 is given by shifting the curve for n = 2 by tt/L, 
and only the solution with n = 1 shows an anomalous 
feature. The curve with n = 1 approaches a point on 
the fc^-axis. This behavior can be checked by setting 
9(k) = tt + k y /k x in Eq. (fT4| with n = 1. It is appro- 
priate to refer to this state with k x = —1/L and k y = 
as the critical state [l5| because this state is neither the 
standing wave nor the edge state The fact that no 
curve crosses the Dirac singularity is consistent to the ab- 
sence of a nontrivial Berry's phase for the standing wave. 
In other words, an electron can not approach the Dirac 
singularity due to the presence of the zigzag edge. As 
a result, an energy gap appears in the spectrum of the 
standing wave. The minimum energy gap is determined 
by the critical state as E gap = 2hvp/L where vf is the 
Fermi velocity. In contrast to the boundary condition for 
the zigzag edge, the boundary condition for the armchair 
edge does not forbid an electronic state cross the Dirac 
singularity point, and therefore the electron can pick up 
a nontrivial Berry's phase. 

Here, we consider pseudospin of the standing wave in 
order to examine the local density of states (LDOS) near 
the zigzag edge. The pseudospin for an eigenstate ^(y) 
is defined by the expectation value of Pauli matrices as 
(<7j) = J <Ji(y)dy (i = x, y, z) where (Ji{y) is a pseudospin 
density defined by (Ji(y) = ^^(y)cr i 'i l (y). The boundary 
condition of Eq. ([7]) means that the pseudospin density is 
polarized into the positive z-axis locally near the zigzag 
edge, that is, a z (0) > and cr x (0) = c y (0) = 0. Ac- 
tually, by putting y ~ into Eq. ||5J), we see that the 
standing wave near the zigzag edge has amplitudes only 
on A-atoms. This polarization of the pseudospin causes 
an anomalous behavior to appear in the LDOS near the 
zigzag edge. To show this, let us first review the LDOS 
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FIG. 3: The curves represent the solutions of parameters k x 
and k y which are allowed by the boundary condition for a 
zigzag nanoribbon. The Dirac singularity is not located on 
the curves, which shows that there is no singularity at the 
parameter space of the standing waves. The inset shows the 
corresponding spectral flow in the energy dispersion relation. 



for a graphene without an edge. Assuming that electrons 
are non-interacting, the bulk LDOS is given by 



P(E) = 



1 



\E\ 



2tt (fo F ) 2 ' 



(15) 



where p(E) is proportional to \E\, which results from the 
Dirac cone spectrum. The LDOS near the zigzag edge 
can be calculated by using the standing wave given in 
Eq. (0). The LDOS has the form, 



Ps(E,y) 



1 



\E\ 



2ir (hv F ) 2 
where R(E, y) is defined as 



R(E,y), 



R(E,y) 



1 



(16) 



(17) 



By performing the integral with respect to the angle 9 in 
Eq. (TP7| . we obtain the analytical result for R(E, y) as 



R(E,y) = 1 



Jo 2 



, \Ey\ 



, 1 2 m\ 

Tiv-p J 



(18) 



where J v (x) is a Bessel function of order v. The LDOS 
is symmetric with respect to E = 0. The bulk LDOS 
in Eq. (jT5j) can be reproduced by setting \y\ —> oo in 
Eq. ([IH| since J„(oo) = 0. In Fig. gj we plot the LDOS 
at \y\ = 0, 1, 2, 3, and 6[nm]. The slope of the LDOS 
depends on the distance from the zigzag edge. At the 
edge, that is, at \y\ — 0, since Jo(0) = 1 and J2(0) = 0, we 
see that R{E, 0) = 1/2, and the slope of p s (E, 0) is half 
of that of p s (E, oo) [= p{E)\. This results from the fact 
that the amplitudes at one of two sublattice disappears 
near the zigzag edge (due to the pseudospin polarization) , 



FIG. 4: The position dependence of the LDOS structure near 
the zigzag edge. The peak structure at E — is due to the 
edge states. 



and only the half of the amplitude in the unit cell can 
contribute to the LDOS. 

In Fig.[U in addition to the LDOS due to the standing 
wave, we p lot the LDOS (at \y\ = 2[nm]) due to the edge 
states The edge states produce a peak at E = 

in the LDOS. The LDOS of the edge states is calculated 
as follows. The wave function of the edge state is given 



by [12J 



K.k 



(r) 



L x 



(19) 



Note that the edge states appear for k x < and y < 
(see Fig. [3]). Since the energy eigenvalue of the edge state 
vanishes, the LDOS for the edge states is written as 



Pe(E,y) 



£ 

fcx<0 



(E - 0) 2 + S 2 ' 



(20) 



where S is a phenomenological parameter representing en- 
ergy uncertainty of the edge states, for which we assume 
6 = 10 mcV. Substituting ^2 k <0 with (L x /2tt) J^^dkx 
in Eq. (|20[) . and performing the integral for k x , we obtain 



Pe{E,y) 



1 



26 



1 



2tt 2 E 2 + S 2 Ay 2 



(21) 



This result has been used in Fig. @]to plot the LDOS of 
the edge states. Note that p e (E,y) decreases as ~ y~ 2 , 
which is a slowly decreasing function compared with the 
exponential decaying wave function of the edge state. It 
is interesting to note that p e (E,y) ~ y~ 2 can be derived 
from the fact that J^ c {p s (E, y) + p e (E, y)}dE does not 
depend on y for a large value of E c . This is a sum rule 
for the total LDOS consisting of the edge states and the 
standing waves, in which the density above (below) the 
Fermi energy at E = must be position independent. 
This argument also shows that the singular behavior of 



5 



the LDOS for the edge states at y = is a consequence 
of the wave function of the standing wave. 

By examining the standing wave near the graphene 
edge, we have seen that the armchair edge preserves a 
non-trivial Berry's phase as in the bulk of graphene, on 
the other hand, such a non-trivial phase is absent for the 
standing wave near the zigzag edge. Here, we consider 
the indication of this result with respect to the trans- 
port in graphene nanoribbons. In the argument for the 
absence of the backward scattering given in the intro- 
duction (see Fig. [T]), we find, by replacing <£> k with the 
standing wave "J/k, that the backscattering amplitude, 
^bs = [1 + U(2w)] \l/-k, is enhanced for zigzag nanorib- 
bons (due to U(2ir) = 1), while it still vanishes for arm- 
chair nanoribbons (due to U(2ir) = — 1). Thus, we can 
expect that the conductance of a zigzag nanoribbon is 
smaller than that of an armchair nanoribbon, if we as- 
sume only scatters with long-range potential. The simu- 
lations performed by Areshkin et al. [18j and Yamamoto 
et al. [19( show indeed that long-range potentials are 
ineffective in causing backscattering in a perfect arm- 
chair nanoribbon. For narrow metallic armchair nanorib- 
bons, the linear energy dispersion originating from the 
Dirac point provides a nearly perfect conduction, and this 
nearly perfect conduction might relate to the nontrivial 
Berry's phase since the linear energy dispersion picks up 
the Dirac singularity. Note, however, that the standing 
wave for the armchair edge is a short-length intervalley 
mode, that is, the electronic wave length of the stand- 
ing wave is order of the lattice constant. As a result, the 
wave function can be very sensitive to a short-range scat- 
tering potential, such as irregular edges [Hj and lattice 
vacancies [2(| • An armchair nanoribbon has both advan- 
tages (non-trivial Berry's phase) and disadvantages (in- 
tervalley mode) for the transport as well as that a zigzag 
nanoribbon has both advantages (intravalley mode) and 
disadvantages (trivial Berry's phase). Therefore, in a 
realistic situation, it is reasonable to consider that a 
nanoribbon does not exhibit an electronic conduction 
comparable to a metallic carbon nanotube. 

It seems that recent tight-binding numerical studies 
show that the conductance of zigzag nanoribbons is more 
robust than that of armchair nanoribbons against edge 
disorders. In the numerical studies, the trans- 

port given by the electrons which are located very close 
to the Fermi level is concerned. For zigzag nanoribbons, 
the propagating modes in each valley contain a single 
one-way excess channel (the states on the curve of n = 1 
which are close to the critical state shown in Fig. [3J . This 
feature causes a perfectly conducting channel to appear 
in the disordered system for the case that impurities do 
not give rise to intervalley scattering [23, 24 1. This fea- 
ture also seems to cause a robust conducting channel to 



companies the corresponding final state with — k x . This 
condition is not satisfied for the single one-way excess 
channel, for which our discussion for the Berry's phase 
is useless. The single one-way excess channel plays a 
dominant role in zigzag nanoribbons of the widths up to 
several nanometers for a limited energy window. (l8j Our 
result for the Berry's phase may be useful in discussing 
transport of zigzag nanoribbons of the widths more than 
10 nm, for which we have several number of channels 
having the final state with —k x . 

Our result may have a relationship to the formation of 
a transport gap observed for graphene nanoribbons. 0^3 
Note that the standing wave near a realistic rough edge 
can be described as the sum of the standing waves for 
zigzag and armchair edges. ll| Since the standing wave 



appear in the edge disordered system. [18|, |20j-|22j For our 



discussion of the backscattering in zigzag nanoribbons, 
we assume that the initial state with momentum k v ac- 



for zigzag edge is a long-length intravalley mode and that 
for armchair edge is a short-length intervalley one, their 
characters do not interfere with each other and can co- 
exist in the standing wave near a rough edge. This is 
also indicative of that a small portion of zigzag edge in a 
rough edge eventually governs the long-length transport 
behavior of a nanoribbon, as is numerically simulated 
by several authors. 2|| 26 1 Indeed, Gallagher et al. 0| 
and Han et al. Q observed that the source-drain gap 
has a strong dependence on ribbon's length. It is also 
important to recognize that there are two distinct local- 
ized states; the edge states and localized states which 
consist of the standing waves. Since the localization 
length of the former is on the order of lattice constant, 
it is difficult to consider that the transport behavior of 
nanoribbons of the widths larger than 10 nm is domi- 
nated by the edge states. [13] In contrast, the latter is 
caused by the scatters, such as a potential disorder cre- 
ated by charge impurities, giving rise to backward scat- 
tering. Since the latter consists of the standing waves, 
there is a possibility that the localization length is the 
order of the width of a ribbon. Han et al. 0] measured 
electron transport in lithographically fabricated nanorib- 
bons of widths 20 < W < 120 nm with rough edge on 
the order of nanometers, and confirmed that the size of 
a transport gap is inversely proportional to the ribbon 
width. 

In conclusion, we have shown that the Berry's phase 
for the standing wave near the zigzag edge is trivial. The 
momentum parameter space for the standing wave near 
the zigzag edge does not include the Dirac singularity, 
which is essential to the absence of a non-trivial Berry's 
phase. As a result, the absence of backward scattering 
mechanism that works well in a metallic carbon nanotube 
can not be used for a zigzag nanoribbon. The absence 
of the Dirac singularity in the momentum space of the 
standing wave results in the peak in the LDOS due to the 
edge states. An observation of the LDOS peak near the 
zigzag edge is a direct evidence of the absence of the Dirac 
singularity in the parameter space of the standing wave. 
In contrast, a non-trivial Berry's phase survives the re- 



6 



flection by the armchair edge. For the standing wave near 
the armchair edge, we obtain the new absence of back- 
ward scattering mechanism in which a real time-reversal 
pair of backward scattered waves cancels with each other. 
This absence of the backward scattering mechanism for 
the armchair edge is not robust against a short-range 
scattering potential since the standing wave is a short- 
length intervalley mode. An armchair nanoribbon has 
both advantage and disadvantage for the transport as 
well as a zigzag nanoribbon. 
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